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Abstract 

Production of D* ± (2010) mesons in diffractive deep inelastic scattering has been 
measured with the ZEUS detector at HERA using an integrated luminosity of 
82 pb _1 . Diffractive events were identified by the presence of a large rapidity 
gap in the final state. Differential cross sections have been measured in the 
kinematic region 1.5 < Q 2 < 200 GeV 2 , 0.02 < y < 0.7, x P < 0.035, /3 < 
0.8, pTiD*^) > 1.5 GeV and \t](D* ± )\ < 1.5. The measured cross sections are 
compared to theoretical predictions. The results are presented in terms of the 
open-charm contribution to the diffractive proton structure function. The data 
demonstrate a strong sensitivity to the diffractive parton densities. 
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1 Introduction 



In ep deep inelastic scattering (DIS) at HERA, final-state hadrons are dominantly pro- 
duced by interactions between virtual photons and incoming protons. Diffractive interac- 
tions, characterized by a large rapidity gap in the distribution of the final-state hadrons, 
have been observed and extensively studied at HERA [11121121111130131311] • The measure- 
ments of the diffractive DIS cross sections [21E1I11E1IZ1 have been quantified in terms of a 
diffractive structure function, F®, defined in analogy with the proton structure function, 
F 2 . The diffractive parton densities, determined from these measurements, are domi- 
nated by gluons. The diffractive process at HERA has often been considered to proceed 
through the exchange of an object carrying the quantum numbers of the vacuum, called 
the Pomeron (IP). In the resolved- Pomeron model JU]; the exchanged Pomeron acts as a 
source of partons, one of which interacts with the virtual photon. In an alternative view, 
the diffractive process at HERA can be described by the dissociation of the virtual photon 
into a qq or qqg state which interacts with the proton by the exchange of two gluons or, 
more generally, a gluon ladder with the quantum numbers of the vacuum fUESEI] • 

Charm production in diffractive DIS, which has also been measured by the HI and 
ZEUS collaborations [THUS], allows quantitative tests of the models due to the sensi- 
tivity of charm production to gluon-initiated processes ^Hj- Calculations based on a 
gluon-dominated resolved Pomeron predict a large charm rate in diffractive DIS [T7lll8j . 
In the two-gluon-exchange models, the rate from the qqg state is similar to that predicted 
by the resolved-Pomeron model, while the rate from the qq state is lower. 

In this analysis, charm production, tagged using D*^ mesons, is studied in diffractive 
interactions identified by the presence of a large rapidity gap between the proton at high 
rapidities and the centrally-produced hadronic system. The luminosity for the present 
measurement is about two times larger than in the previous ZEUS study [T^- The increase 
in luminosity and an improved rapidity acceptance in the proton direction allow a more 
detailed comparison with the model predictions in a wider kinematic range. The open- 
charm contribution to the diffractive proton structure function is measured for the first 
time. 



2 Experimental set-up 

The analysis was performed with data taken from 1998 to 2000, when HERA collided 
electrons or positrons 1 with energy E e = 27.5 GeV with protons of energy E p = 920 GeV 

1 Hereafter, both e + and e~ are referred to as electrons, unless explicitly stated otherwise. 
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yielding a centre-of-mass energy of 318 GeV. The results are based on the sum of the 
e~p and e + p samples, corresponding to integrated luminosities of 16.4 ± 0.3 pb -1 and 
65.3 ± 1.5 pb" 1 , respectively. 

A detailed description of the ZEUS detector can be found elsewhere [T^]. A brief out- 
line of the components that are most relevant for this analysis is given below. Charged 
particles are tracked in the central tracking detector (CTD) [20], which operates in a 
magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD consists 
of 72 cylindrical drift chamber layers, organised in nine superlayers covering the polar- 
angle 2 region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks 
is a{p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [2I| consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/ \f~E for electrons 
and a(E)/E = 0.35/\/E for hadrons, with E in GeV. The timing resolution of the CAL 
is better than 1 ns for energy deposits greater than 4.5 GeV. 

In 1998-2000, the forward plug calorimeter (FPC) [22] was installed in the 20 x 20 cm 2 
beam hole of the FCAL, with a small hole of radius 3.15 cm in the centre to accom- 
modate the beam pipe. The FPC increased the forward calorimetric coverage by about 
1 unit of pseudorapidity to t]<5. The FPC consisted of a lead-scintillator sandwich 
calorimeter divided longitudinally into electromagnetic and hadronic sections that were 
read out separately by wavelength-shifting fibers and photomultipliers. The energy res- 
olution, as measured under test-beam conditions, was a(E)/E = 0Al/y/E © 0.062 and 
a(E)/E = 0.65/ VE © 0.06 for electrons and pions, respectively, with E in GeV. 

The position of electrons scattered at a small angle with respect to the electron beam 
direction was measured using the small-angle rear tracking detector (SRTD) [221 • The 
luminosity was determined from the rate of the bremsstrahlung process ep —>■ ejp, where 
the photon was measured with a lead-scintillator calorimeter [21] located at Z = —107 m. 



2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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3 Kinematics and reconstruction of variables 



The four-momenta k, k' and P label the incoming electron, outgoing electron and the 
incoming proton, respectively, in DIS events: 

e(k) + p(P) — > e(k') + anything. 

To describe the kinematics of DIS events, any two of the following invariants can be used: 

r\2 2 ( i i/\2 Q 2 P-q 2 Q 2 (l — x) 

Q =~q = ~(k - k ) ; x = ——; y = -— -; W = - 



2P-q y P-h' x 

where Q 2 is the negative square of the four-momentum q carried by the virtual photon, x 
is the Bjorken scaling variable, y is the fraction of the electron energy transferred to the 
proton in its rest frame, and W is the centre-of-mass energy of the photon-proton system. 
The scattered electron was identified using an algorithm based on a neural network |25j . 
The hadronic final state was reconstructed using combinations of CTD tracks and energy 
clusters measured in the CAL and FPC to form energy-flow objects (EFOs) |51l7|l26j . The 
kinematic variables were reconstructed using the double-angle method |2*7j . 

To describe the diffractive process ep — ► eXp, where X is the hadronic final state origi- 
nating from the dissociation of the virtual photon, two additional variables were used: 

* x P = Q 2 +w% ) where Mx is the invariant mass of the system X. This variable is the 
fraction of the incoming proton momentum carried by the diffractive exchange; 

* @ = x~p~ = q 2 +m' 2 • I n an interpretation in which partonic structure is ascribed to the 
diffractive exchange, j3 is the longitudinal momentum fraction of the exchange that is 
carried by the struck quark. 

The above expressions neglect the proton mass. The square of the four-momentum trans- 
fer at the proton vertex, t, was not measured; thus all results are implicitly integrated 
over this variable, which was assumed to be zero in the expressions for x P and f3. 

The mass of the diffractive system X was calculated from EFOs using: 

mi = (x>v - (e ^,<) 2 - (e p y) 2 - (E a* 

where the sum i runs over the EFOs not associated with the scattered electron. 

The process studied in this paper is ep — > eXp — > e{D* ± X')p, in which the system X 
includes at least one D* ± meson. The latter was reconstructed using the mass- difference 
method [2H1 i n the decay channel D* + —> D ir+ followed by D° — ► K~ir + (+c.c), where 
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ir s indicates the "slow" pion. The fractional momentum of the D*^ 1 in the photon-proton 
system is defined as 



W ' 

where p*(D* ± ) is the D* ± momentum in the photon-proton centre-of-mass frame. 



4 Models of diffractive charm production 

In the resolved-Pomeron model, proposed by Ingelman and Schlein ^U]; the exchanged 
Pomeron is assumed to be a object with a partonic structure. The diffractive cross section 
factorises into a Pomeron flux factor, describing the probability to find a Pomeron in the 
proton; the Pomeron's parton density function (PDF), specifying the probability to find 
a given parton in the Pomeron; and the interaction cross section with the parton. Within 
this model, open charm is produced in diffractive DIS via the boson-gluon-fusion (BGF) 
process, where the virtual photon interacts with a gluon from the Pomeron (Fig.[IJi). The 
HERA measurements of the inclusive diffractive differential cross sections were found to 
be consistent with the resolved-Pomeron model with a Pomeron structure dominated by 
gluons. For x P > 0.01, an additional contribution from Reggeon exchanges, carrying the 
quantum numbers of a p, u, a or / meson, was found to be sizeable jlj. A combined fit 
of the Pomeron parton densities to the HI and ZEUS inclusive diffractive DIS measure- 
ments |HEll2niE] and to the ZEUS data on diffractive dijet photoproduction [20] has been 
made by Alvero et al. (ACTW) [21]. The Pomeron flux factor was assumed to be of the 
Donnachie-Landshoff form [22] arid only data satisfying x P < 0.01 were used. To fit the 
Pomeron parton densities, five functional forms (labelled A, B, C, D and SG) were used. 
It was found that only gluon-dominated fits (B, D and SG) were able to describe both 
the DIS and photoproduction data, while the quark-dominated fits (A and C) underesti- 
mated the photoproduction data significantly. Therefore, only the gluon-dominated fits 
are compared to the data in Section 8. The fit results have been interfaced to the pro- 
gram HVQDIS [22] to calculate cross sections for diffractive charm production in DIS [18J, 
both to leading and next-to-leading order (NLO) in QCD. In this analysis, the ACTW 
NLO predictions were calculated setting the charm-quark mass m c = 1.45 GeV and the 
renormalisation and factorisation scales hr = fip = \/Q 2 + 4rn% as in [IB]. The Peter- 
son fragmentation function (with e = 0.035 [21]) was used for the charm decay. The 
probability for charm to fragment into a D* ± meson was set to f(c — > D* + ) = 0.235 |35j . 

The two- gluon- exchange models consider fluctuations of the virtual photon into qq or 
qqg colour dipoles that interact with the proton via colour- singlet exchange; the simplest 
form of which is a pair of gluons j2E] • The virtual-photon fluctuations into cc (Fig. [IJ)) 
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and ccg states (Fig. Hfc) can lead to diffractive open-charm production. At high x P val- 
ues, quark exchanges are expected to become significant. Thus, the two-gluon-exchange 
calculations are expected to be valid only at low x P values (x P < 0.01). In recent calcu- 
lations |HZllI21ISHllSn] , the cross section for two-gluon exchange is related to the square of 
the unintegrated gluon distribution of the proton which depends on the gluon transverse 
momentum, kx, relative to the proton direction. In the "saturation" model jSHHHII, the 
calculation of the qqg cross section is performed under the assumption of strong kx order- 
ing of the final-state partons, which corresponds to kjfi *C k^\ The parameters of the 
model were tuned to describe the total photon-proton cross section measured at HERA. 
Alternatively, in the model of Bartels et al. [23II21EH1 , configurations without strong kr 
ordering are included in the qqg cross-section calculation and the minimum value for the 
final-state-gluon transverse momentum, fc^ u *, is a free parameter. The sum of the cc and 
ccg contributions in the saturation model and the model of Bartels et al. are hereafter re- 
ferred to as SATRAP and BJLW, respectively. Both the SATRAP and BJLW predictions 
were calculated using the MC generator RAPGAP 2.08/06 jlj, the proton PDF param- 

eterisation GRV94HO [12], m c = 1.45 GeV and = fiF = \Jp1,t + 4m^, where 
p Ct T is the transverse momentum of the charm quark. Such scale form was used because 
RAPGAP does not provide the form used for the ACTW predictions. The probability for 
open charm to fragment into a D*^ meson was set to f(c — > D* + ) = 0.235. In the BJLW 
calculation of the ccg component, the value of the parameter /c™* was set to 1.5 GeV |43| . 



5 Acceptance calculation 

To study trigger and selection efficiencies, two MC programs, RAPGAP and RIDI 2.0 |44j . 
were used to model the final states in the process ep — *> eXp — > e{D* ± X')p. 

The RAPGAP generator was used in the resolved-Pomeron mode, in which charm quarks 
are produced via the leading-order BGF process of Fig.^i. The higher-order QCD correc- 
tions were simulated using the colour-dipole model implemented in ARIADNE 4.03 |45| . 
The LUND string model jlH] as implemented in JETSET 7.4 j37j was used for hadroni- 
sation. The charm-quark mass was set to the default value of 1.5 GeV. The diffractive 
sample was generated assuming a gluon-dominated Pomeron, with a parameterisation 
from the HI Collaboration called "HI fit 2" [T3|. The Reggeon (meson) component of 
the parameterisation was not used. 

The RIDI generator is based on the two-gluon-exchange model developed by Ryskin |44j . 
To simulate the gluon momentum density, the GRV94HO proton PDF parameterisation 
was used. Final-state parton showers and hadronisation were simulated using JETSET 
and the charm-quark mass was set to the default value of 1.35 GeV. First-order radiative 
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corrections were included in the simulation although their effects were negligible. The 
cc and ccg components were generated separately and later combined in the proportion 
16% : 84% which provided the best description of the ft distribution of the data. 

The RAPGAP MC sample was used to evaluate the acceptance. Three MC samples 
were used to estimate the model dependence of the acceptance corrections: the RIDI 
MC sample, a sample generated with RAPGAP using parton showers as implemented in 
LEPTO 6.1 |48J to simulate the higher-order QCD corrections, and a sample generated 
with RAPGAP using the Pomeron PDF parameterisation "HI fit 3" [T3j . 

To estimate the non-diffractive DIS background and to measure the ratio of diffractive 
to inclusive D* ± production (see Section l8~2*j) . two MC generators were used: RAPGAP 
in the non-diffractive mode for the nominal calculations and HERWIG 6.301 jlH] as a 
systematic check. The RAPGAP parameters used were the same as those used in the 
ZEUS measurement of the inclusive DIS -D* ± cross sections [SU]. To generate charm pro- 
duction via the leading-order BGF process with HERWIG, the CTEQ5L proton PDF 
parameterisation and m c = 1.5 GeV were used. Hadronisation in HERWIG is simulated 
with a cluster algorithm |52j . 

In this analysis, the final-state proton was not detected. To estimate and subtract the 
contribution from the diffractive processes where the proton dissociates into a system N, 
ep — > eXN — > e{D* ± X')N ', four MC generators were used: DIFFVM [53] for the nominal 
calculations and RAPGAP, PHOJET [S3] and EPSOFT 2.0 [SS] for systematic checks. 
The DIFFVM MC program provides a detailed description of the proton-dissociative fi- 
nal state. The mass spectrum, Mjv, of the system N is generated as a superposition of 
N* + resonances and a continuum having the form da/dM N oc M^ 2 e \ The default 
parameter value e = 0.0808 [SH] was used. In the RAPGAP simulation of proton dis- 
sociation, the proton splits into a quark and di-quark and the Pomeron is assumed to 
couple only to the single quark. The spectrum follows a 1/Mjv distribution. In PHO- 
JET, Mm is calculated from the triple-Pomeron kinematics [S3] and an approximation of 
the low-mass-resonance structure. In EPSOFT, the M^-spectrum generation relies on a 
parameterisation of the pp —* pN data. 

The generated events were passed through the GEANT-based [HZ] simulation of the ZEUS 
detector and trigger. They were reconstructed by the same program chain as the data. 
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6 Event selection and D reconstruction 

6.1 Trigger and DIS selection 

Events were selected online with a three-level trigger fHlEB]- At the first level, events 
with an electron candidate in the EMC sections of RCAL or BCAL were selected jHH]- I n 
the latter case, a coincidence with a track originating at the nominal interaction point was 
required. At the second level, the non-ep background was further reduced by removing 
events with CAL timing inconsistent with an ep interaction. At the third level, events 
were fully reconstructed and selected by requiring a coincidence of a scattered-electron 
candidate found within the CAL and a D*^ candidate reconstructed in the nominal de- 
cay mode using charged tracks measured by the CTD. The requirements were similar 
to, but looser than, the offline cuts described below. The efficiency of the online D*^ 
reconstruction, determined relative to an inclusive DIS trigger, was above 95%. 

The following criteria were applied offline to select DIS events: 

• an electron with energy above 10 GeV; 

• the impact point of the scattered electron on the RCAL lies outside the region 26 x 14 
cm 2 centred on the beamline; 

• 40 < 8 < 65 GeV, where 5 = JZiC^i — Pz,%) and the sum runs over the EFOs from the 
hadronic system and the energy deposited by the identified electron; 

• a vertex position |Z vtx | < 50 cm. 

The events were restricted to the kinematic region 1.5 < Q 2 < 200 GeV 2 and 0.02 < y < 
0.7. 

6.2 D* ± reconstruction 

Charged tracks with p T > 0.12 GeV and \rj\ < 1.75 were selected. Only tracks assigned 
to the primary event vertex and with hits in at least three superlayers of the CTD were 
considered. Two oppositely charged tracks, each with p? > 0.5 GeV, were combined to 
form a D° candidate. The tracks were alternately assigned the mass of a charged kaon 
and a charged pion and the invariant mass of the track pair, M(Kir), was calculated. 
Only D° candidates that satisfy 1.81 < M(Kir) < 1.92 GeV were kept. Any additional 
track, with p? > 0.12 GeV and charge opposite to that of the kaon track, was assigned the 
pion mass and combined with the D° candidate to form a D*^ candidate with invariant 
mass M{Ktvtx s ). The D*± candidates were required to have pt{D*) > 1.5 GeV and 
\rj(D*)\ < 1.5. 
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In the distribution of the mass difference, AM = M{Kt[tt s ) — M(Kn), for selected D*^ 
candidates, a clear signal at the nominal value of M(D* ± ) — M(D°) was observed (not 
shown). The combinatorial background under this signal was estimated from the mass- 
difference distribution for wrong-charge combinations, in which both tracks forming the 
D° candidates have the same charge and the third track has the opposite charge. The 
number of reconstructed D* ± mesons was determined by subtracting the wrong-charge 
AM distribution after normalising it to the AM distribution of D*^ candidates with the 
appropriate charges in the range 0.15 < AM < 0.17 GeV. The subtraction, performed 
in the range 0.1435 < AM < 0.1475 GeV, yielded an inclusive signal of 4976 ± 103 D*± 
mesons. 

6.3 Selection of diffractive events 

Diffractive events are characterised by the presence of a large rapidity gap between the 
proton at high rapidities and the centrally-produced hadronic system. To select such 
events, the following two requirements were applied: 

• -Efpc < 1-5 GeV, where -Efpc is the energy deposited in the FPC; 

• ^max < 3, where ?7 max is the pseudorapidity of the most-forward EFO measured without 
using FPC information and with energy above 400 MeV. 

This selection is illustrated in Fig. where the distribution of r/ max is shown for D*^ 
mesons obtained after the wrong-charge-background subtraction. The data are compared 
to the ?7 max distributions of the non-diffractive RAPGAP and HERWIG MC samples 
and to the sum of the non-diffractive and diffractive RAPGAP MC. In Fig. |2k, the 
distributions are shown for events with any -Eppc value. The large peak at ?? max ~ 3.5 
corresponds to non-diffractive events in which the proton remnant deposits energy around 
the beam direction. On the low side of the peak, the contribution from non-diffractive 
interactions exhibits an exponential fall-off, leaving an excess at low values of 77 max which 
is populated predominantly by diffractive events. Figure shows that the requirement 
Efpc < 1-5 GeV strongly suppresses the contribution from non-diffractive interactions. 
Requiring r/ max < 3 in addition reduces the remaining non-diffractive background and 
ensures a gap of at least two units of pseudorapidity with respect to the edge of the 
forward calorimetric coverage (see Section |2J). 

The selected events were analysed in terms of the diffractive variables x P , j3 and Mx- To 
account for the restriction imposed by the ?7 max < 3 requirement, a cut of x P < 0.035 was 
applied. In addition, a cut of (3 < 0.8 was also used because diffractive charm production 
in DIS is strongly suppressed at large (3 values due the dominant contribution of events 
with small Q 2 and large Mx values. 
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Figure El shows the AM distribution after the above cuts. The number of D*^ 1 after the 
wrong-charge-background subtraction is 253 ±21. 

Figure |U shows the number of reconstructed D*^ mesons in bins of the variables pr(D* ± ), 
r](D* ± ), x(D* ± ), p, x P , \og(Mx), log(Q 2 ) and W. The data are compared to the diffractive 
RAPGAP and RIDI simulations (normalised to the data). Both simulations reproduce 
the shapes of the data. 

6.4 Subtraction of the proton-dissociative contribution 

Diffractive events with proton dissociation can pass the -Efpc < 1-5 GeV and ?7 max < 3 
requirements if the major part of the proton-dissociative system escapes undetected down 
the forward beampipe. The proton-dissociative contribution was determined from the 
distribution of E FP q for events selected with relaxed reconstruction cuts and with- 
out cutting on -Efpc- To ensure a gap of at least two units of pseudorapidity between 
the proton-dissociative system, tagged by the FPC, and the system X, a requirement of 
Vmax < 1-75 was applied. Figure El compares the -Eppc distribution for these events to the 
distributions of the diffractive RAPGAP and proton-dissociative DIFFVM MC samples. 
The MC samples were combined in the proportion providing the best description of the 
£"fpc distribution, and their sum was normalised to the data. Using the normalisation 
factors obtained for the two MC samples, the proton-dissociative contribution was calcu- 
lated for the nominal diffractive selection described in Section l6~3l The proton-dissociative 
contribution was determined to be 16% with negligible statistical uncertainty; the sys- 
tematic uncertainty was obtained as follows, where the effects of each source are shown 
in parentheses: 

• the parameter b, regulating the shape of the Mjy continuum distribution in the DIF- 
FVM MC simulation, was varied between 0.7 and 1.5 (^39%); 

• uncertainties in the low-mass resonance structure and other details of the simulation 
of the proton-dissociative system were estimated by using the PHOJET, RAPGAP 
and EPSOFT MC generators 

• a shift of ±10% due to the FPC energy-scale uncertainty (1qj%); 

• a larger area, including the FPC and neighbouring FCAL towers, was used to tag the 
proton-dissociative system (—2.7%). This check is sensitive to the high-M^v proton- 
dissociative contribution and to details of the FPC and FCAL simulation. 

These systematic uncertainties were added in quadrature separately for the positive and 
negative variations to determine the overall systematic uncertainty of ±4.1%. The proton- 
dissociative contribution of (16 ± 4)% was assumed to be independent of all kinematic 
variables and was subtracted from all measured cross sections. 
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7 Systematic uncertainties 



The systematic uncertainties of the measured cross sections were determined by changing 
the selection cuts or the analysis procedure in turn and repeating the extraction of the 
cross sections [HD]- The major sources of the systematic uncertainty were as the follows, 
where effects on the integrated cross section are shown in parentheses: 

• the selection of inclusive DIS events (133%). Variations were made in the cut on the 
scattered-electron energy, the RCAL box cut, the 5 cut and the vertex-position cut. 
In addition, both Q 2 and y were determined using the eS method [HI] rather than 
using the DA method; 

• the selection of D*^ candidates and background estimation (^37%). The minimum 
transverse momentum for the K and ir candidates was raised and lowered by 25 MeV. 
For the slow pion, tt s , the minimum transverse momentum was raised and lowered by 
10 MeV. The signal region for M(D°) was loosened to 1.80 < M(D°) < 1.93 GeV and 
that of the AM distribution was widened to 0.143 < AM < 0.148 GeV. The AM 
background-normalisation region was varied by 5 MeV; 

• the selection of diffractive events The requirements on ?7 max and Eppc were 
varied by ±0.2 units and ±0.5 GeV, respectively; 

• a shift of ±3% due to the CAL energy-scale uncertainty 

• a shift of ±10% due to the FPC energy-scale uncertainty (^3%); 

• the model dependence of the non- diffractive contribution (—6.6%). This uncertainty 
was estimated using the HERWIG sample; 

• the model dependence of the acceptance corrections (1 7 ; 4 %). This uncertainty was es- 
timated using the RIDI MC sample, the RAPGAP sample generated with the LEPTO 
parton showers and the RAPGAP sample generated with the "HI fit 3" parameteri- 
sation of the Pomeron structure function. 

These systematic uncertainties were added in quadrature separately for the positive and 
negative variations to determine the overall systematic uncertainty of -fi,2%- These esti- 
mates were also made in each bin in which the differential cross sections were measured. 

The normalisation uncertainties in the luminosity measurement (±2.2%) and the D*^ and 
D° branching ratios (±2.5% [H2]) were not included in the systematic uncertainty. The 
uncertainty arising from the subtraction of the proton-dissociative background, quoted 
separately, is ±4.1%/0.84 = ±4.9%. 
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8 Results 



8.1 Cross sections 

The differential D* ± cross sections for any given variable £ were determined using: 

da _ N(D*) (1 - f pd ) 
d£ A C B A£ ' 

where N(D*) is the number of D*^ mesons in a bin of width A£, A is the acceptance 
for that bin, £ is the integrated luminosity, B is the product of the D* + — ► D°ir+ and 
.D — > K~ir + branching ratios (0.0257 [62J), and (0.16) is the fraction of the proton- 
dissociative background discussed in Section lfi~4l 

Using the overall acceptance of 19.4%, the cross section for diffractive D*^ production 
in the kinematic region 1.5 < Q 2 < 200 GeV 2 , 0.02 < y < 0.7, x P < 0.035, (3 < 0.8, 
Pt{D* ± ) > 1.5 GeV and |r 7 (Z)* ± ) | < 1.5 is 

<? ep ^ e D*±x'p — 521 ± 43(stat.)l5g(syst.)±26(p.diss.) pb, 

where the last uncertainty arises from the subtraction of the proton-dissociative back- 
ground 3 . 

In the case of Reggeon exchanges, open charm can be produced in the BGF process if 
the exchanged-meson PDF contains gluons. The Reggeon contribution to diffractive D*^ 1 
production in the measured kinematic range was estimated to be less than 6% using 
RAPGAP with the Pomeron and meson PDF parameterisations "HI fit 2" or "HI fit 
3". The contribution is less than 0.5% for x P < 0.01; it increases with x P , contributing 
about 12% in the last bin. The Reggeon contribution, which is smaller than the statistical 
uncertainty of the measurement, was neglected. 

Figure El (Table [TJ shows the differential cross section as a function of x P . The data are 
compared with the ACTW NLO predictions, calculated with the gluon-dominated fit B, 
the SATRAP predictions and the BJLW predictions. All three models agree with the data 
within experimental uncertainties below x P = 0.01. For larger x P values, the ACTW and 
SATRAP models agree with the data whereas the BJLW prediction underestimates the 
measured cross sections as expected (see Section HJ. 

The differential cross sections as functions of priD^), ?7(D* ± ), log(M|), x^D^), /3, 
log(/3), log(<5 2 ) and W were measured for x P < 0.01 and x P < 0.035 (Tables 121 and |HJ). 
Figure [7| compares the differential cross sections measured for x P < 0.01 with the ACTW, 

3 The diffractive Z3* ± cross section was also calculated in the kinematic regions in which previous 
measurements |f 4lll5j were reported and was found to be consistent. 
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SATRAP and BJLW predictions. In Figs. Hand H the ACTW and SATRAP predictions 
are compared with the differential cross sections measured for x P < 0.035. 

The two-gluon-exchange BJLW model predictions, obtained with the cutoff value /c™* = 
1.5 GeV tuned using the HI measurement [Tl], describe the differential cross sections in 
the range x P < 0.01 both in shape and normalisation. Using the value = 1.0 GeV 
(2.0 GeV), the model predictions significantly overestimate (underestimate) the data in 
this range (not shown). 

The two-gluon-exchange saturation model (SATRAP) predictions reproduce the shapes 
and the normalisations of the differential cross sections measured in both x P ranges. 

The ACTW NLO predictions, obtained with the gluon-dominated fit B, describe the data 
reasonably well in both x P ranges. Using other gluon-dominated fits, the predictions 
significantly overestimate (fit D) or underestimate (fit SG) the data (not shown). The 
quark-dominated fits A and C were excluded by the previous ZEUS measurement [To] . 

8.2 Ratio of diffractive to inclusive J D* ± production 

The ratio of diffractively produced D mesons to inclusive D mesons, Rd, was mea- 
sured for x < 0.028. This limit is the product of the x P and (5 requirements imposed for 
the diffractive D* ± sample. The ratio of diffractive to inclusive DIS D*^ production is 
then defined by 

_ a ep ^ eD *±x'p(x P < 0.035, (3 < 0.8) 
o ep ^ eD *± Y {x < 0.028) 

Sources of systematic uncertainty in the ratio measurement were studied in a similar 
manner to those for the cross-section measurements. There is a cancellation between the 
common systematic uncertainties originating from the selection of inclusive DIS events, 
the selection of D*^ candidates and the background estimation. An additional contri- 
bution originates from the model dependence of the acceptance corrections used in the 
evaluation of the inclusive DIS D*^ cross sections. This systematic uncertainty was esti- 
mated using the inclusive RAPGAP MC sample generated with LEPTO parton showers 
instead of the ARIADNE higher-order QCD corrections and with the HERWIG MC sam- 
ple. 

The ratio measured in the kinematic region 1.5 < Q 2 < 200 GeV 2 , 0.02 < y < 0.7, 
PriD**) > 1.5 GeV, |?7(D* ± )| < 1.5 and x < 0.028 is 

R D = 6.4±0.5(stat.)^(syst.)+^(p.diss.) %. 

The value is consistent with previous measurements performed in similar kinematic ranges 

CSJ. 
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Figure ITU1 f Table E|) shows the ratio measured as a function of pxiD*^), ^(D* ± ), x(-D* ± ), 
log(Q 2 ) and W. The measured Rd shows no dependence on Q 2 , W or x^D^). The 
relative diffractive contribution is larger at small pr(-D* ± ) and in the backward direction 
(negative rj(D )). The NLO QCD predictions for the ratio of diffractive to inclusive DIS 
D* ± production were obtained using ACTW NLO fit B for the diffractive predictions and 
the HVQDIS program with the CTEQ5F3 [SI] proton PDF for the inclusive predictions. 
Parameters in both calculations were set to the values discussed in Section 0] The NLO 
QCD predictions reproduce the measured Rd values and the trends observed for the Rd 
distributions measured as functions of pt^D*^) and ?7(D* ± ). 



8.3 Open-charm contribution to the diffractive proton structure 
function F2 

Neglecting contributions from Z-boson exchange and the longitudinal structure function, 
the open-charm contribution to the diffractive structure function of the proton can be 
related to the cross section, measured in the full D*^ kinematic region, by 

2f(c^D*+) dx P d(3dQ 2 ~ Q*P {1 V+ 2>* 2 > x *)' W 

In order to estimate F®^ ,cc ^ the differential cross section was measured as a function of 
\og{(3) for different regions of Q 2 and x P (Table EJ). Extrapolation factors of the measured 
cross sections to the full pt{D*^) and ^(D^) phase space were estimated using the ACTW 
NLO fit B predictions. The factors were about five for low-x P bins and two for high-x^. 
bins. 

In each bin, J?®^ ,cc W as determined using the formula 



!,meas 

■pD(S),cc, R n 2 \ _ a ep-+eD*±X>p p D{3),cc ( a n 2 

r 2 measl^) ; X JP,«J — t,ACTW 2 ACTWlW' Vi ) x p,«;i 

<J ep->eD* ± X'p 



where the cross sections a 1 in bin % are those for pr(D* ± ) > 1.5 GeV and |?7(D* ± )| < 1.5. 
The functional form of -^actw: ca l cma t e d using Eq. (JTJ), was used to quote the results 
for i^ 2 D ( 3 ) ,cc a t convenient values of Q 2 and x P) i close to the centre-of-gravity of the 
bin. 

The measured F^ 3 '' cc values are listed in Table El with their experimental uncertainties. 
Using ACTW NLO fit D had no significant effect on the measured values. Other sources 
of extrapolation uncertainties are small compared to the experimental uncertainties |5T)] . 

Figure E] shows the quantity function of log(/5) for different Q 2 and x P 

values. In all cases, XpF®^^ rises as (3 decreases. The curves show the theoretical 



13 



XjpF 2 obtained using the ACTW NLO calculations with fit B, D and SG. The fit 
B prediction generally agrees with the data. The fit D (SG) prediction overestimates 
(underestimates) the measured x P F®^' cc at low (3. 

9 Summary 

Diffractive D*^ production has been measured in the kinematic region 1.5 < Q 2 < 
200 GeV 2 , 0.02 < y < 0.7, x P < 0.035, (3 < 0.8, prC **) > 1-5 GeV and M-D**)] < 1.5. 
The cross section integrated over this kinematic region is 521 ± 43(stat.)+^(syst.) ± 
26(p.diss.) pb. Differential cross sections have been compared to the predictions of dif- 
ferent diffractive models. The ACTW NLO predictions, based on parton densities of the 
Pomeron obtained from combined fits to the inclusive diffractive DIS and diffractive di- 
jet photoproduction measurements at HERA, describe the results reasonably well in the 
whole ge if the gluon-dominated fit B is used. The predictions of the two-gluon- 

exchange saturation model also reproduce the shapes and normalisations of the differential 
cross sections in the whole x P range. The predictions of the two-gluon-exchange BJLW 
model describe the cross sections measured for x P < 0.01, if a minimum value for the 
transverse momentum of the final-state gluon of kj^ = 1.5 GeV is used. 

The ratio of diffractive D* ± production to inclusive DIS D*^ production has 
been measured to be Rd = 6.4 ± 0.5(stat.)^Q;7(syst.)^o;3(p.diss.) %. The ratio Rd shows 
no dependence on W, Q 2 or x{D* ± ). The relative contribution from diffraction is larger 
at small Pt(D ) and in the backward direction (negative ?7(D* ± )). The NLO QCD 
predictions reproduce the measured Rd- 

The open-charm contribution, p®^ ,cc ^ to the diffractive proton structure function has 
been extracted. For all values of Q 2 and x P , F% 00 rises as (3 decreases. The results have 
been compared with the theoretical F®^ ,cc obtained using the ACTW NLO calculations 
with the gluon-dominated fits B, D and SG. The data exclude the fits D and SG, and 
are consistent with fit B. This demonstrates that the data have a strong sensitivity to 
the diffractive parton densities, and that diffractive PDFs in NLO QCD are able to 
consistently describe both inclusive diffractive DIS and diffractive charm production in 
DIS. 
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x P bin 


da/dx P (nb) 


, 0.003 


28.0 ±4.9^2 


0.003 , 0.006 


25.4±4.7^i 


0.006 , 0.010 


18.6 ±3.6^| 


0.010 , 0.020 


13.7 ±2.2^;| 


0.020 , 0.035 


13.7±2.4l|° 



Table 1: Differential cross section for diffractive D production as a function 
of Xp . The first and second uncertainties represent statistical and systematic un- 
certainties, respectively. The overall normalisation uncertainties arising from the 
luminosity measurement (±2.2%), from the D* ± and D° branching ratios (±2.5%) 
and from the proton- dissociative background subtraction (±4.9%) are not indicated. 
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PtGD**) bin ( GeV) 


da/dp T ( J D* ± ) (pb/GeV) 




x P < 0.01 


x P < 0.035 


1.5 , 2.4 


161 ± 29tl 8 9 


307 ± 50+^ 


2.4 , 3.3 


66 ± 111? 


151 ± 20±i| 


3.3 , 4.2 


19 ± 5+ 2 


70 ± 11±| 


4.2 , 5.4 


10 ± 3±{ 


26 ± 5±1 


5.4 , 10.0 




2.8 ± 0.9±g;i 


T]{D*^) bin 


da/dr/( J D* ± ) (pb) 




x P < 0.01 


x r < 0.035 


-1.5 , -0.9 


124 ± 26+}^ 


212 ± 36±|? 


-0.9 , -0.3 


104 ± 19+* 4 


213 ± 31+^ 


-0.3 , 0.3 


78 ± tftl 1 


195 ± 29±p 


0.3 , 0.9 


37 ± 13±f 2 


125 ± 28±|| 


0.9 , 1.5 


55 ± 20±l\ 


134 ± 36±H 


log(Mi/GeV 2 ) bin 


d ( j/dlog(M|-/GeV 2 ) (pb) 




x P < 0.01 


x P < 0.035 


1.00 , 1.44 


89 ± 21±2i 


94 ± 23+ 22 


1.44 , 1.88 


195 ± 35±H 


201 ± 38±| 


1.88 , 2.32 


200 ± 29± 2 , 4 


382 ± 


2.32 , 2.76 


A7±25t\l 


284 ± 54+^ 


2.76 , 3.20 




286 ± 65i}^ 


x(D* ± ) bin 


da/dx( J D* ± ) (pb) 




x F < 0.01 


x P < 0.035 


0.00 , 0.16 


185 ± Glttl 


429 ± 1071^ 


0.16 , 0.32 


252 ± 76t 7 5 2 


788 ± 1351^ 


0.32 , 0.48 


446 ± 85j^ 


864 ± 1341^ 


0.48 , 0.64 


376 ± 75t% 


726±119i^ 


0.64 , 1.00 


92 ± 21+^ 8 


221 ± 38+ 2 ^ 



Table 2: Differential cross sections for diffractive D production as a function 
of px{D* ± ), viD*^), log(Mjf) and x{D* ± ). The first and second uncertainties 
are statistical and systematic, respectively. The overall normalisation uncertainties 
arising from the luminosity measurement (±2.2%), from the and D° branching 
ratios (±2.5%) and from the proton-dissociative background subtraction (±4.9%o) 
are not indicated. 
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(3 bin 


dtr/d/3 (pb) 




x P < 0.01 


x P < 0.035 


0.00 , 0.10 


1252 ± 2031}™ 


4153 ± 410+ 2 ^ 


0.10 , 0.20 


419 ± 94j1j2 


654 ± I25tnt 


0.20 , 0.30 


244 ± 54+^ 


311 ± 69+H 


0.30 , 0.45 


100 ± 35+^ 


91 ± 39lg 


0.45 , 0.80 


27 ± H+5 4 


33 ± 13±^ 5 


\og((3) bin 


da/dlog(/3) (nb) 




x P < 0.01 


x P < 0.035 


-3.0 , -2.0 




115 ±33 + ;g 


-2.0 , -1.5 


105 ± 28l| 


392 ± 58l^ 


-1.5 , -1.0 


124 ± 25^ 


272 ± Alttj 


-1.0 , -0.5 


141 ± 22l} 2 


203 ± 28+26 


-0.5 , -0.1 


65 ± 16+£ 


56 ± 18l^ 7 


log(Q 2 / GeV 2 ) bin 


da/dlog(Q 2 /GeV 2 ) (pb) 




x P < 0.01 


x P < 0.035 


0.17, 0.6 


276 ± 51±H 


534 ± 871^ 


0.6 , 1.0 


140 ± 29tf 5 


324 ±5111 


1.0 , 1.3 


106 ± 27l* 


342 ± 50lg 


1.3 , 1.55 


103 ± 25l}8 


225 ± 43l^ 


1.55 , 2.3 


17 ± 7i£ 


41 ± 13+J 6 


W bin ( GeV) 


da/dW (pb/ GeV) 




x P < 0.01 


x P < 0.035 


50 , 92 


0.45 ±0.1418;" 


1.53±0.35+°; 2 j 


92 , 134 


1.48 ± 0.29±g;lf 


3.36 ± 0.49±^;|f 


134 , 176 


1.63 ±0.29^i? 


3.68 ± 0.491^50 


176 , 218 


1.25±0.29+°; 2 ^ 


2.43 ± 0.44l8;37 


218 , 260 


0.50 ± 0.3318; 22 


0.95 ± 0.48l8;ts 



Table 3: Differential cross sections for diffractive D production as a function 
of (3, \og{(3), log(Q 2 ) and W . The first and second uncertainties are statistical and 
systematic, respectively. The overall normalisation uncertainties arising from the 
luminosity measurement (±2.2%), from the D*+ and D° branching ratios (±2.5%) 
and from the proton-dissociative background subtraction (±4.9%o) are not indicated. 
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PrCD**) bin ( GeV) 


Rd (%) 


1.5 , 2.4 


8.5 ±1.5123 


2.4 , 3.3 


6.3 ±0.9i£ 2 . 


3.3 , 4.2 


5.5±0.9i°J 


4.2 , 5.4 


4.3 ± 0.91^ 


5.4 , 10.0 


2.5 ± 0.8l°; 2 


rjiD*^) bin 


Rd (%) 


-1.5 , -0.9 


11.2 ±2.0^1 


-0.9 , -0.3 


8.6 ± l-3±?;^ 


-0.3 , 0.3 


6.8±l.lt°j 


0.3 , 0.9 


4.4 ± 1.0±g;t 


0.9 , 1.5 


4.4 ± 1.2±2;? 


x(D* ± ) bin 


Rd (%) 


0.00 , 0.16 


5.0±1.3i^ 


0.16 , 0.32 


6.2 d= i.i±i:o 


0.32 , 0.48 


6.4 il.Olg;} 


0.48 , 0.64 


7.4±1.2i?;jj 


0.64 , 1.00 


9.6 ± 1.7±§;| 


log(Q 2 / GeV 2 ) bin 


Rd (%) 


0.17 , 0.60 


7.9±l.3i2;? 


0.60 , 1.00 


5.8±0.9i?;^ 


1.00 , 1.30 


8.1 ± 1.2±g;| 


1.30 , 1.55 


7.8 ± 1.6l°; 2 


1.55 , 2.30 


3.6 ± 1.2±g;t 


W bin ( GeV) 


(%) 


50 , 92 


5.1 ± 1.2±2;| 


92 , 134 


6.6 ± 1.0l?;o 


134 , 176 


7.7±l.li°| 


176 , 218 


7.4 ±1.41^ 


218 , 260 


4.4 ±2.3lJ;g 



Table 4: Ratio of diffr actively produced D*^ mesons to inclusive D*^ meson 
production as a function of ' pt{D*^) , ^(D*^), x(D* ± ), log(Q 2 ) and W . The first 
and second uncertainties are statistical and systematic, respectively. The over- 
all normalisation uncertainties arising from the luminosity measurement (±2.2%), 
from the D*^ and D° branching ratios (±2.5%) and from the proton- dissociative 
background subtraction (±4.9%) are not indicated. 
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log(/3) bin 


da/dlog(/3), x P < 0.01 (pb) 




1.5 < Q 2 < 10 GeV 2 


10 < Q 2 < 200 GeV 2 


-2.0 , -1.5 


107 ± 28+| 




-1.5 , -1.0 


114 ±25+3° 




-1.0 , -0.5 


62 ± 16^ 4 


80 ± 15±| 


-0.5 , -0.1 




61 ± 16+] 3 


log(/3) bin 


do-/dlog(/3), 0.01 <x F < 0.035 (pb) 




1.5 < Q 2 < 10 GeV 2 


10 < Q 2 < 200 GeV 2 


-3.0 , -2.0 


96 ± 31+J 




-2.0 , -1.5 


142 ± 431^ 


141 ± 30±^ 


-1.5 , -1.0 




106 ± 25+j 3 , 


-1.0 , -0.5 




52 ± ntll 



Table 5: Differential cross section for diffractive D production as a function 
of log(/5) for different regions of Q 2 and x P . The first and second uncertainties 
are statistical and systematic, respectively. The overall normalisation uncertainties 
arising from the luminosity measurement (±2.2%), from the D*+ and D° branching 
ratios (±2.5%) and from the proton-dissociative background subtraction (±4.9%o) 
are not indicated. 
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F 2 D(3)CC ", x P = 0.004 


(3 


Q 2 = 4 GeV 2 


Q 2 = 25 GeV 2 


0.020 


1.34 ± 0.35j£|t 




0.050 


0.92 ± 0.20+°:?3 




0.200 


0.20 ± 0.05^ 


2.14 ±0.40lgijj 


0.500 




0.89 ± 0.23±g:ii 


Ff (3)c5 , x F = 0.02 


/? 


Q 2 = 4 GeV 2 


Q 2 = 25 GeV 2 


0.005 


0.20 ± 0.071°;^ 




0.020 


0.17 ±0.05^ 


1.87 ±0A0toil 


0.050 




0.50 ±0.121^9 


0.200 




0.18 ±0.06^ 



Table 6: The measured charm contribution to the diffractive structure function 
of the proton, F®^' 00 , f or different values of (3, Q 2 and x F . The first and second 
uncertainties are statistical and systematic, respectively. The overall normalisation 
uncertainties arising from the luminosity measurement (±2.2%), from the -D* ± and 
D° branching ratios (±2.5%) and from the proton-dissociative background subtrac- 
tion (±4.9%o) are not indicated. 
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(b) (c) 

Figure 1: Modelling charm production in diffractive ep scattering: (a) boson-gluon 
fusion in the resolved-Pomeron model, (b) cc and (c) ccg states in the two-gluon- 
exchange model. 
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ZEUS 




Figure 2: Numbers of reconstructed D*^ mesons (dots) as a function of 77 max 
for DIS events with (a) any -E^pc values and (b) -Efpc < 1-5 GeV . The solid 
histogram shows the sum of the non-diffractive RAPGAP MC (hatched area) and 
the diffractive RAPGAP MC. The sum was normalised to have the same area as 
the data. The dashed histogram shows the non-diffractive HERWIG MC. 
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ZEUS 



• ZEUS 98-00 
x Ip < 0.035, P < 0.8 
— Wrong-charge background 




M(KtoO - M(Kjc) (GeV) 



Figure 3: The distribution of the mass difference, AM = M(K7nr s )— M(Ktt), for 
candidates (dots) in events with r] ma _ x < 3, E FPC < 1-5 GeV, x P < 0.035 and 
(3 < 0.8. The histogram shows the AM distribution for wrong-charge combinations. 
Only candidates from the shaded band were used for the differential cross- 
section measurements. 
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Figure 4: Numbers of reconstructed D*^ mesons (dots) in bins of pt^D*^), 
TjiD**), x{D*±), f3, x P , log(M|), log(Q 2 ) and W. The RAPGAP (solid histogram) 
and the mixed cc and ccg RIDI (dashed histogram) MC samples, normalized to the 
data, are shown for comparison. 
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Figure 5: The measured energy in the FPC for events with ?7 max < 1.75 (dots). 
The dashed histogram is the single- diffractive RAPGAP MC sample and the dotted 
histogram is the proton- dissociative DIFFVM MC sample. The solid histogram is 
the sum of both diffractive and proton-dissociative MC samples normalised to the 
data. 
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cut 
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Figure 6: Differential cross-section da/dx P for diffractive D*^ production for the 
data (dots) compared with the ACTW NLO (solid histogram), SATRAP (dashed 
histogram) and BJLW (dotted histogram) predictions. The shaded area shows the 
effect of varying the charm-quark mass in the ACTW NLO prediction. The inner 
error bars indicate the statistical uncertainties, while the outer ones correspond to 
statistical and systematic uncertainties added in quadrature. The overall normal- 
isation uncertainties arising from the luminosity measurement (±2.2%), from the 
D*^ and D° branching ratios (±2.5%o) and from the proton-dissociative background 
subtraction (±4.9%) are not indicated. 
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Figure 7: Differential cross sections for diffractive D*^ production with x P < 0.01 
for the data (dots) compared with the ACTW NLO (solid histogram), SATRAP 
(dashed histogram) and BJLW (dotted histogram) predictions. The shaded area 
shows the effect of varying the charm quark-mass in the ACTW NLO prediction. 
The cross sections are shown as a function o/pt(-D* ± ) j ^(-D* 1 * 1 ); log(M|-) and f3. 
The inner error bars indicate the statistical uncertainties, while the outer ones cor- 
respond to statistical and systematic uncertainties added in quadrature. The over- 
all normalisation uncertainties arising from the luminosity measurement (±2.2%), 
from the -D* ± and D° branching ratios (±2.5%) and from the proton-dissociative 
background subtraction (±4.9%o) are not indicated. 
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Figure 8: Differential cross sections for diffractive D production with x P < 
0.035 for the data (dots) compared with the ACTW NLO (solid histogram) and 
SATRAP (dashed histogram) predictions. The shaded area shows the effect of vary- 
ing the charm-quark mass in the ACTW NLO prediction. The cross sections are 
shown as a function o/pt(-D* ± ) j 1 '?(-D* ± ); log(^x) an ^ x{D* ± ). The inner er- 
ror bars indicate the statistical uncertainties, while the outer ones correspond to 
statistical and systematic uncertainties added in quadrature. The overall normal- 
isation uncertainties arising from the luminosity measurement (±2.2%), from the 
D*^ and D° branching ratios (±2.5%) and from the proton-dissociative background 
subtraction (±4.9%o) are not indicated. 
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Figure 9: Differential cross sections for diffractive D production with x P < 
0.035 for the data (dots) compared with the ACTW NLO (solid histogram) and 
SATRAP (dashed histogram) predictions. The shaded area shows the effect of vary- 
ing the charm-quark mass in the ACTW NLO prediction. The cross sections are 
shown as a function of (3, log(/3) ; log(Q 2 ) and W. The inner error bars indicate 
the statistical uncertainties, while the outer ones correspond to statistical and sys- 
tematic uncertainties added in quadrature. The overall normalisation uncertainties 
arising from the luminosity measurement (±2.2%), from the D*^ and D° branching 
ratios (±2.5%) and from the proton-dissociative background subtraction (±A.Q%o) 
are not indicated. 
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Figure 10: The measured ratio of diffractively produced D* ± mesons to inclu- 
sive D*^ meson production (dots). The ratio is shown as a function o/pr(D* ± ) J 
•qip^), x{D* ± ), log(Q 2 ) and W . The inner error bars indicate the statistical 
uncertainties, while the outer ones correspond to statistical and systematic uncer- 
tainties added in quadrature. The histogram corresponds to the NLO QCD predic- 
tion where the shaded area shows the effect of varying the charm-quark mass. The 
overall normalisation uncertainty arising from the proton- dissociative background 
subtraction (±.A.Q%) is not indicated. 
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Figure 11: The measured charm contribution to the diffractive structure func- 
tion of the proton multiplied by function of (3 for different 
values of Q 2 and x P (dots). The inner error bars indicate the statistical uncer- 
tainties, while the outer ones correspond to statistical and systematic uncertainties 
added in quadrature. The overall normalisation uncertainties arising from the lu- 
minosity measurement (±2.2%), from the D*^ and D° branching ratios (±2.5%) 
and from the proton-dissociative background subtraction (±4.9%) are not indicated. 
The curves correspond to the ACTW model prediction; the shaded area shows the 
effect of varying the charm-quark mass. 
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